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We compute the finite temperature correction to the induced fermion number for fermions coupled 
to a static SU(2) chiral background, using the derivative expansion technique. At zero temperature 
the induced fermion number is topological, being the winding number of the chiral background. 
At finite temperature however, higher order terms in the derivative expansion give nontopological 
corrections to the winding number. We use this result to show that the standard cancellation of the 
fractional parts of the fermion number inside and outside the bag in an SU(2) x SU(2) hybrid chiral 
bag model of the nucleon does not occur at nonzero temperature. 
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I. INTRODUCTION 



When second quantized fermions interact with a classical background field with a nontrivial topology (e.g, solitons, 
vortices, monopoles, Skyrmions etc) the resulting quantum states of the system can possess fractional eigenvalues. 
This intriguing phenomenon, known as fermion number fractionization, has applications ranging from particle, nuclear 
and atomic physics to condensed matter systems p], 0, B 0, [| pi |7| §]. At zero temperature, the fermion number 
of the vacuum is a topological quantity (up to spectral flow effects), and is related to the spectral asymmetry of the 
relevant Dirac operator, which counts the difference between the number of positive and negative energy states of the 
fermion spectrum. Rigorous mathematical results, such as index theorems and Levinson's theorem, show that the 
fractional part of the vacuum fermion number is a topological invariant; i.e, it depends only on the asymptotic values 
of the background field, and is invariant under local deformations of the background || [l0| pT| . This topological 
character of the fermion number is important for various applications in model field theories, such as soliton models 
for the nucleon since it allows the fermion number to be kept fixed in a variational calculation that minimizes the 
^> . energy |l3|, However, at finite temperature, the induced fermion number is generically nontopological, and 
moreover is not a sharp observable p"o| fl6| , p"7[ [Tc| . Several explicit examples of finite temperature fermion number in 
kink 15, sigma model 15, j^J and monopole Jl5|, |l) backgrounds have been analyzed in detail. In this paper, 



we compute the finite temperature nontopological corrections to the induced fermion number for a (3+1) dimensional 
<-j ', SU(2) chiral background, and show explicitly that the finite temperature correction is nontopological. 

The model of second quantized fermions chirally coupled to classical, heavy scalar and pseudoscalar fields has 
. £h ', applications in particle and nuclear physics, where it is used in models to describe the low energy spectrum of QCD 
such as in chiral quark soliton models of hadrons jl^, ^2], |23| and in the hybrid chiral bag model of the nucleon 

SHI- At 

zero temperature, the induced fermion number of fermions interacting with a static chiral field coincides 
with the winding number of the field, and this result may be interpreted in terms of a topological current density. 
Physically, the topological zero temperature contribution to the fermion number corresponds to the virtual dipole 
pairs of vacuum polarization, which do not live long enough to thermalize. The nontopological finite temperature 
corrections correspond to the real qq pairs in the thermal plasma which are sensitive to the details of the single particle 
spectrum. 

This paper is organized as follows. In section II, we define the finite T induced fermion number and review why 
it is generically nontopological. In Section III we use the derivative expansion to compute higher order T dependent 
corrections to the winding number of the chiral field. As a nontrivial check we show that it vanishes in the T — > limit. 
We conclude in section IV and discuss the implications of the existence of nontopological T dependent corrections for 
the hybrid chiral bag model for the nucleon. 



II. FINITE TEMPERATURE INDUCED FERMION NUMBER 



The induced fermion number is an expectation value of the second quantized fermion number operator, N = 
\ J dx For a given static classical background field configuration, the second quantized fermion field operator 

^> can be expanded in a complete set of states of the single particle Dirac Hamiltonian H. At zero temperature, the 
fermion number is a vacuum expectation value (N)o = (0|iV|0), and by a standard Fock space calculation H is related 
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to the spectral asymmetry of the Dirac Hamiltonian 



(N)o — — — (spectral asymmetry) 



dEa(E) sign(£) (1) 



Here a(E) is the spectral function, or density of states, of the Dirac Hamiltonian H : 

a{E) = - Im Tr ( 1 . ) (2) 

7r \H — E — te J 

Thus the fermion number of the vacuum (N)o is essentially the difference between the number of positive and 
negative energy states of the Dirac spectrum. When fermions interact with topologically nontrivial backgrounds 
(solitons, vortices, monopolcs, skyrmions), or are constrained by certain boundary conditions (e.g., in the chiral bag), 
the spectral asymmetry can be nonintegral, because of distortions in the Dirac sea [Q, 0, 0|. 
At nonzero temperature, T, the induced fermion number is a thermal expectation value : 



dEa(E) tanh (^-J (3) 

where /3 = i. This finite temperature expression (0) reduces smoothly to the zero temperature expression (1) in 
the zero temperature limit f3 — > oo. The expression ml) also emphasizes the fact that the spectrum of single particle 
states for fermions interacting with the classical background field is completely determined by the single particle 
Dirac Hamiltonian H of the system and has nothing to do with temperature. The temperature only influences the 
probability of each energy level being occupied. The finite temperature induced fermion number (N)t splits naturally 
into a zero temperature piece (N)q, given by (Q) and a temperature dependent correction 



(N) T = {N) + dEa(E)sign(E)n(\E\) (4) 

J — OO 

Here n(E) is the Fermi-Dirac distribution function 

and we have used the simple identity tanh(/3_E/2) = 1 — 2n(E). Thus the temperature enters only in the second 
term of (|j) via the Fermi-Dirac distribution n(E). While {N} is topological, the finite temperature correction to the 
induced fermion number, being dependent on the Fermi weighting factor n(|.E|), is more sensitive to the details of the 
fermion spectrum, which in turn depends on the precise shape of the background field. Therefore, the generic situation 
is that the finite temperature correction is nontopological. However, in certain special cases, when the background is 
such that the Dirac Hamiltonian has a quantum mechanical SUSY, then the odd part of the spectral function a(E) 
[which is what is needed in (^)] is itself topological, and therefore the finite temperature correction is also topological 
even though it is temperature dependent. The simplest cases where this happens are for a kink background in 1+1 
dimensions |]l5[ |l9| , or for a monopole background!^, |2l| . 

The physical interpretation of the separation in (Q)~is the following Jig] . The zero temperature fermion number (N)o 
represents the vacuum polarization of the fermions in the presence of the background; it is temperature independent 
because the virtual dipolcs do not live long enough to come to thermal equilibrium. The temperature dependent 
correction corresponds to the plasma response of the system in the presence of the background and is thus sensitive to 
the details of the spectrum. The finite temperature induced fermion number (and the rms fluctuations of the fermion 
number) have been computed for sigma model backgrounds in (1+1) and (2+1) dimensions (l^, |l6[ [Tt], ^8). In this 
paper we consider second quantized fermions chirally coupled to an SU(2) background in (3+1) dimensions. 
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III. INDUCED FERMION NUMBER IN A CHIRAL FIELD 

A. The Model 

Consider an isodoublet of fermions interacting in an SU(2)^ x SU(2)# chiral invariant way with a classical, static 
chiral background in (3+1) dimensions 

C = i§7 M 9 M * - * (er + i7 5 f • if) * 

= i^j^d^ - exp (vy^r ■ (j^j * (6) 

Here a is a scalar field and if is a triplet of pseudoscalar (pion) fields. The Pauli matrices, r, are the generators of 
SU(2), and m is the mass scale in the theory. The Goldstone boson fields (a, if) are constrainted to lie on the "chiral 
circle", er 2 + if • if — m , and thus they can be parametrized by a unitary matrix U as 



U 



exp (if -(f)] = — (er + if ■ if) (7) 
V / m 



Thus 



d = mcos(|0|) 

if = m sm(\(j)\)if (8) 

where the hat denotes a unit vector in isospin space. The model (^) is used in various quark-soliton models of baryons 
wherein the baryon is described in terms of quarks interacting with a soliton of the classical chiral field |l^, ^] . The 
model (||) also finds applications in extensions of the MIT bag model of the nucleon, to incorporate chiral symmetry; 
in these "hybrid" chiral bag models B |l3), quarks confined inside a boundary are chirally coupled to scalar and 
pseudoscalar fields (mesons) at the bag boundary with appropriate boundary conditions to maintain continuity of the 
axial vector current. 

At zero temperature, the induced fermion number for the model (Q) can be computed using the derivative expansion 
technique and is given by ^, |22|, pi[ 

Wo = ^j J ' #xJ ik Tr{d i UU*d j Urfd k Urt) (9) 

Thus, the zero temperature induced fermion number is simply the winding number of the group clement t/(x). 
The winding number is a topological invariant, determined by the asymptotic behavior of U(x). Invoking Lorentz 
invariance, the topological charge density (the integrand of (N)o) can be interpreted as the zeroth component of a 
current density 

B„ - ^ e^Tr (d v Utfd p Utfd°Utf) (10) 

This current is often called a topological current because it is conserved, c^B^ 1 = 0, without use of the equations 
of motion. In fact, the current density B^ in ([l0]) has corrections involving higher order derivatives of the field 
U. However, these higher derivative corrections are total derivatives and thus do not contribute to an integrated 
quantity such as the fermion number (N)o (H[ ^2|. At nonzero temperature, experience with similar models in lower 
dimensions |l6| , p"8| suggests that the fermion number density Bo receives higher order derivative corrections that are 
not total derivatives and that are nontopological. This expectation will be confirmed by an explicit calculation in 
the next section for the 3 + 1 dimensional chiral quark model (|). This appearance of nontopological corrections at 
finite temperature is completely consistent with the fact that the Dirac Hamiltonian for this system has no quantum 
mechanical SUSY. Indeed, numerical results [[12], |25| for the spectrum of fermions in the background of an SU(2) chiral 
field U, in the "hedgehog" ansatz, show that the quark energy levels are highly sensitive to the shape (in particular, 
the length scale) of the hedgehog profile function. By the arguments of the previous section, this means that the finite 
temperature corrections in (|j) are necessarily nontopological. 

B. The Derivative Expansion Calculation 

To calculate the induced fermion number at nonzero temperature we adopt the general methodology used in p5[ . 
Using the ie prescription in the spectral density (|^) we rewrite (N)t as contour integral using the Dirac resolvent 
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Here z is an arbitrary complex number, and H is the single particle Dirac Hamiltonian corresponding to (g|) in the 
fermion sector. Thus 

where C is the contour (—00 + ie, +00 + ie) and (+00 — ie, —00 — ie) in the complex energy plane. For a Hermitian 
Hamiltonian H , the poles and branch cuts of the resolvent lie on the real axis. On the other hand, the Matsubara poles 
z n = (2n + l)«7rT of the tanh(/3z/2) function lie on the imaginary axis. Thus we have two equivalent representations 
for (N)t ■ deforming the contour C around the poles and cuts of R(z) leads to an integral representation, and 
deforming C around the discrete z n leads to a summation representation. The integral representation is just the 
familiar Sommerfeld- Watson transform of the summation representation. Thus the computational strategy is: 

(i) Obtain an expression for the Dirac resolvent R(z). In most cases this requires an approximation since the Dirac 
equation in the presence of complicated backgrounds can rarely be solved exactly, especially in higher dimensions. In 
this paper we use the derivative expansion to obtain an approximate expression for R{z). 

(ii) Insert the resulting expression for R{z) into ([12]) to obtain an expression for (N)t either as an integral, or as 
a sum over the Matsubara frequencies z n . 

We stress that these two steps are independent. In particular, the derivative expansion calculation has nothing to 
do with the temperature - it is simply a means to obtain approximate information about the spectral properties of 
the Hamiltonian in the presence of the static b ackg round. 

In the derivative expansion approximation || one assumes that derivatives of the background field (here the 

chiral field U) are much smaller than the mass scale to in the theory, 

|W|<m (13) 

Since we are considering a static background configuration all these derivatives are spatial derivatives. So the back- 
ground field is assumed to be approximately constant on the length scale of the fermion Compton wavelength. First, 
note that to compute (N)t we need the odd part of the spectral density u{E), or equivalently the even part (in z) of 
the resolvent R(z). The even part of the resolvent (O) can be written as 



Tr 



1 



H-z 



'■(//tTT 1 -:) CL4) 



even 



Our strategy is to develop a systematic derivative expansion for this even part of the resolvent. From the Lagrangian 
@, the corresponding Dirac Hamiltonian for the fermions in the chiral background U is 

H = -il ® 7°7 • V + mcos(\(j)\)I ® 7 + imr ■ tt ® 7°7 5 sin(|$) 

= -il (g> a 1 (g> a 1 d t + y [(U + U^)(g> a 3 ® I + i{U - <8> a 2 ® i] (15) 

where we have used the definitions (0) and (|) and I is the 2x2 identity matrix. The first term is the kinetic energy of 
the fermions, the second is the chiral invariant fermion-meson interaction. We choose the gamma matrix convention 
as in ||: 7 = er 3 ® I, j l — ia 2 ® a 1 , j 5 = a 1 % I. The square of the Hamiltonian (|l5| ) is 

H 2 = I® I® I (-V 2 + m 2 ) + y [i (diU - diU*) ®a 3 ®a l - (diU + d. t U r ) ® a 2 ® a 1 ] 

= H 2 + V (16) 
where the interaction term V is 

V = 7? [* (9iU - 8^) ®a 3 ® a 1 - (d.U + d^) ® a 2 ® a 1 } (17) 

and Hq — (— V 2 + to 2 ) is the square of the free Dirac Hamiltonian. Since V involves one factor of m and one derivative 
of U, it is small compared to the natural scale to 2 of Hq, owing to the derivative expansion condition (|l3|). Hence the 
even part of the resolvent (n3) can be systematically expanded in powers of V: 



Tr 



1 



H- 

where the free propagator A is 



Tr{H AVAyAVA) + Tt(HAVAVAVAVA) - Tr(i?AFAFAFAFAFA) + • • • (18) 



A = (19) 

— V 2 + TO 2 — Z 2 
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The trace in (|18|) involves a matrix trace over the Dirac indices, a trace over the internal isospin indices, and a 
functional trace. The leading term in this expansion involves three factors of V, since the Dirac and Isospin traces 
involving H and either one or two factors of V give zero. Note that in computing these traces the operators H and A 
act on everything to the right of them. We thus need to know the action of A on the perturbation V(~x). We make 
use use of the following operator identity: 



1 



-V(x) - V(x) 



1 



-V 2 + a 
+4 {didjV^didj 



2 (W(x) 
1 



1 



(-V 2 + a 2 ) 2 



v 2 v( x ; 



i 



(-V 2 + a 2 ) 



2VS 



v (v 2 v(x; 



(-V 2 

1 



a 2 ) 2 



(-V 2 +a 2 ) 3 



(20) 



where the parentheses around the derivatives on the RHS of ( |20| ) indicate that the gradient operator acts on V(x) 
only. At any given order of the derivative expansion we collect together the required number of derivatives of the field 
U (all possible combination of derivatives for a given order), as given by repeated application of the above identity, 
noting that V(x) itself already involves one derivative of U. 

The leading term in the derivative expansion involves three factors of V from the first term in ( |l8| ) . This corresponds 
to all factors of V pulled out, with A having no action on them (corresponding to the first term in (po|)): 



Tr 



1 



H 



(3) 



= -Tt{HAVAVAVA) {3) 

= rr 2 ^ / d 3 x^ k ^{ud^d 3 ud k u^ 

lD7r(m z — z z ) 2 J 



(21) 



Substituting this for (N)t in (|12|), we perform the contour integral as a sum over Matsubara modes, and find the 
finite T contribution at this order of the derivative expansion to be 



n— — oo 




m T 



16tt [to 2 + ((2n+ 1)ttT) 2 ] 5 

3tTTO 4 T \ 

2 [to 2 + ((2n+ 1)ttT) 2 ]^ / 



d 3 x e ijk Tr (Ud t U^d j Ud k U^) 



(22) 



Notice that this leading-order term is simply the zero temperature fermion number (N)q, multiplied by a tempera- 
ture dependent prefactor. In the low temperature limit, this leading-order term ([22]) reduces smoothly to the zero 
temperature fermion number (iV)o in (|^): 



Here we have used the general low temperature expansion 

1 m 1 - 2p T(p- |) (2mT) 1 -P 



T E 



n— — oo 



(((2n + 1)ttT) 2 + to 2 )p 20? r(p) mr(p) 



e t + 



(23) 



(24) 



Thus, in the zero temperature limit the leading derivative expansion term involving three derivative of U gives the 
entire zero temperature answer (Bh, with temperature dependent corrections that vanish exponentially as T — > 0. 

The leading-order expression ( p2[) is still topological since it is simply the zero temperature answer multiplied by 
a smooth function of temperature. The nontopological contributions to (N)o come from the next higher term in the 
derivative expansion which involves five derivatives of U (the fourth order terms in the derivative expansion give zero 
because the functional trace over the propagators turn out to be zero). At fifth order, each of the first three terms 
in the derivative expansion ( |l8| ) contributes. Using the operator identity (|o|) we collect together terms involving five 
derivatives of U from each of the three terms in (|l8|). For the last term involving five factors of V(x) this corresponds 
to taking all the V's out with the A's having no effect on them. For the first two terms we need to go to higher 
derivative terms in ( p0| ) to collect the five derivative terms. After a considerable amount of manipulation we have 
found the following compact expression for this fifth order contribution: 



Tr 



H-z 



l(5) 



to 6 + 6m 4 z 2 
2567t(to 2 — z 1 



d 3 x e ijk Tr 



6 



-4 (udiU^djUd k U^(WU • W f )} 



+ - 



Jjk 



Tr 



2567t(to 2 — z 2 ) 2 
-4<9 4 {[/^[^(VdfcE/ • W f ) - Ud.jU\VU ■ Vd, 

-~6d t iu ] djU{WU^ ■ Vd k U) - U^d 3 U{Vd k U^ ■ VU)\ 

-10 d l {udjU^dkUV 2 ^ - tfdjUdkU^u} 

In obtaining ( p5|) we have used the following identity 



5V • {diU^djUdkUtVU - fyV djtf d k UVU^\ 



(25) 



Tr 



^UditfdiUdkU* - vWdiUditfdkU 



= -Tr 



UditfdjUdkU^VU* ■ W) - U ] d.Ud^ ~d k U{VU • W f ) 



which can be easily proved using the unitarity of U, the cyclicity of the trace and the fact that ■ VU = 

(Vcr • Vcr + V7r a • V7r a )7, where U is parametrized as in (0). 

From (^5|) we see that the fifth-order term in the even part of the resolvent contains two types of contributions: 
first, those that cannot be expressed as a total derivative (the first two terms), and second those which are total 
derivatives. The total derivative terms go to zero after doing the spatial integrals. The remaining terms in (|2^) are 
clearly nontopological and cannot be expressed in terms of the algebraic structure of the winding number density in 
(^) . These terms give a nonzero contribution to the spatial integral and depend sensitively on the profile of the field 
U, not just on the asymptotic behavior of U. Substituting the fifth-order resolvent expression (p5) into the general 
expression ([^J) for the finite temperature fermion number, it is a simple matter to evaluate the z integral. This gives 
the fifth order contribution to the finite T induced fermion number: 



(N)P 



= ~T E 



to 6 - 6m 4 ((2n + 1)ttT) 2 
256tt [to 2 + ((2n+ 1)ttT) 2 ]^ 



d 3 xe ijk Tr UdiUdjU^VU ■ Vd k U^) ~ d l U' l d j U(WU^ ■ Vd k U)\ 
-A^Ud^djUdkU^VU ■ W f )}] (26) 



Putting everything together the finite temperature induced fermion number {N)x of fermions coupled to a chiral field, 
up to fifth order in the derivative expansion, is the sum of (p2|) and (|2q): 



(N) T = 



T 



OO 

E 



37TTO 4 T 



21 1 



„ = _ oc 2[m 2 + ((2n+l)^T) 

to 6 - 6m 4 ((2n + 1)ttT) 2 



24tt 2 



d 3 x e tjk Tr (UditfdjUdktf) 



E 



-oc 



256tt [m 2 + ((2n + 1)ttT) 



d 3 cce ij ' fe Tr 



-■1- 



■ VdkU^) - d^djUiVU^ ■ Vd k U)\ 
{udiU^djUdkU^WU- W f )} +••■ (27) 



This is the main result of this paper. The first term is proportional to the (topological) winding number and 
as shown in (|22| J23| ) , the temperature dependent prefactor smoothly reduces to one (plus exponentially decaying 
corrections) as T — > 0, and we recover the well-known zero temperature answer that the fermion number is equal to 

the winding number of the background field U. The second set of terms in ( p7| ) (coming from (N)j! ), are clearly 
nontopological. A nontrivial consistency check is that this nontopological contribution must vanish in the T — > limit , 
since the T = fermion number is topological and comes entirely from the third-order derivative expansion term. 
This is shown as follows: using the low temperature expansion (|24|), the T dependent prefactor of the nontopological 
piece in (|27|) behaves as (keeping leading and next-to-leading terms) 



T 



OO 

E 

i— — OO 



to 6 - 6to 4 ((2?i + 1)ttT) 2 
2567r[TO 2 + ((2n + l)7rT) 2 ]i 



256tt 



/2toT / 



V2 V ^F T(9/2) mT(9/2) 

a / m- 6 T(3) (2mT)- 5 / 2 /T 
V20Fr(7/2) mT(7/2) 



3840ttT 4 



T 

1 - 6— 

to 



(28) 
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The important thing to note here is that the leading constant pieces cancel, leaving the leading contribution as an 
exponentially suppressed term which vanishes as T — ► 0. This should be contrasted with the low temperature limit of 
the prefactor of the third order term in (pz]), where a constant term survives, yielding the windin g nu mber, as in ([23]). 
The cancellation in (|2q ) of the leading constant pieces in the prefactor of the fifth order term in (|27|) means that this 
fifth order contribution (which is nontopological) to the fermion number vanishes at T = 0. This is a stringent check 
on the fact that the zero temperature fermion number is known to be topological. We stress that this cancellation is 
highly nontrivial because the compact fifth order expression (^6|) comes about by combining three different terms in 
the fifth order of the derivative expansion ( |l8| ) . 



IV. CONCLUSIONS AND DISCUSSION 



We have presented an explicit derivative expansion computation of the finite temperature corrections to the induced 
fermion number for fermions coupled to a static SU(2) chiral background. This calculation illustrates the splitting 
of the induced fermion number into a zero temperature piece, which is topological (here it is the winding number 
of the chiral field), and a finite temperature correction which is nontopological. That the finite T induced fermion 
number will have a nontopological contribution was argued on general grounds in jl?], [l8| where sigma models in 
(1+1) and (2+1) dimensions were analyzed. It was shown in that generically there exist nontopological contributions 
to the induced fermion number, except for very special background fields for which the Dirac Hamiltonian has a 
quantum mechanical SUSY relating the positive and negative energy spectra. The Dirac spectrum in the presence 



of a "hedgehog" chiral background has been calculated numerically in |12, and clearly shows the absence of any 
quantum mechanical SUSY, as well as the sensitive dependence of the spectrum on the scale of the background. 
Since the fermion spectrum is not symmetric and is sensitive to the scale of the background (i.e., the length scale of 
the hedgehog profile function), the finite T corrections must also be sensitive to this scale because of the presence 
of the Fermi weighting factor n{E) in (|J). Thus the finite temperature corrections must be nontopological in this 
model. We have given an analytic calculation that confirms this numerical expectation and moreover yields the precise 
compact form M27J) of the next order in the derivative expansion for the fermion number at finite temperature. Our 
calculation of the finite T induced fermion number was done in the derivative expansion limit which assumes that the 
spatial derivatives of the background field are much smaller than the fermion mass scale in the theory. A nontrivial 
consistency check on our result is that the nontopological finite T contribution, coming from the fifth order in the 
derivative expansion, has a prefactor that vanishes exponentially in the T — > limit. 

An interesting implication of our result is for the SU^(2) x SU,r(2) hybrid chiral bag model of the nucleon. This 
model consists of quarks confined in a three dimensional spherical bag, of radius R, with appropriate boundary 
conditions so that the quark current vanishes outside the bag, ensuring confinement |}], ^3|. Inside the bag the quarks 
are free and at the bag surface they obey the chiral boundary condition 

[ij ■ n + exp (if • n 7 5 6>)] *] bag boundary = (29) 

where \& is the quark field and h is the unit normal vector at the bag boundary. From the point of view of inside 
the bag, the parameter 9 is simply a boundary condition parameter. The boundary condition ( |29| ) is consistent with 
chiral symmetry in the sense that it ensures that the axial current is continuous across the bag surface. Goldstone 
and Jaffe computed the zero temperature fermion number with N c = 3 quarks inside the bag |s| 



B in (0) = 1 - - 

7T 



- sin 29 
2 



(30) 



The first term is simply the contribution of the valence quarks and the fractional piece comes from the contribution of 
the polarized Dirac sea. This polarization arises because the quark spectrum differs from the free spectrum owing to 
the nontrivial boundary condition (|2^). The topological nature of this zero temperature fermion number is illustrated 
by the fact that ( |30| ) does not depend on the bag radius R. 

In the hybrid chiral quark bag model || [u|, on the outside of the bag the quarks are coupled to meson (Skyrme) 
fields, which are usually taken to be spherically symmetric, with the pseudoscalar pion fields in (j^) being in the 
"hedgehog ansatz", <f> = f9(r). Outside the bag the quarks are coupled to a chiral field exactly as in the Lagrangian 
(0). So the induced fermion number comming from the region outside the bag can be calculated by computing the 
space integral in (0), using the hedgehog ansatz for U, from R to infinity H. This gives 



Bout — — 
7T 



0(R) - -sin20(R) 



(31) 
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The topological nature of this contribution to the zero temperature fermion number is illustrated by the fact that 
( |3l| ) only depends on the boundary values of the background field, not on its detailed shape. Next, the inside and 
outside regions are joined by identifying the boundary conditionparameter 9 with the value of the hedgehog profile 
field at the bag boundary: 9 — 9(R). Then combining (^o|) and (|ol| ) we observe that the fractional parts of ( p0|) and 
( |3l| ) cancel, leaving the total fermion number of this hybrid chiral bag equal to one, the baryon number H. One can 
thus model the nucleon as a chiral bag defect in a Skyrme background. Clearly, this cancellation depends crucially 
on the topological nature of the induced fermion number, both inside and outside the bag. 

We now ask whether this cancellation persists at nonzero temperature. From our results in this paper it is clear that 
the contribution from outside the bag, where the quarks are coupled to a chiral field, will be nontopological at nonzero 
temperature. The fifth-order term in (p7|), when integrated over the external region from R to spatial infinity, clearly 
yields a temperature dependent contribution that is highly sensitive to the details of the profile function 9(r) appearing 
in the hedgehog field U, and not just to the boundary value 9{R). On the inside of the bag, the finite temperature 
induced fermion number is also nontopological. This follows from the fact that the detailed quark spectrum, with 
the boundary condition (p9|), is sensitive to the bag radius R. This can be deduced from the explicit solutions in a 
grand spin basis (see jtl| for a phase shift analysis of the related external problem with boundary condition (p9|)). 
Now suppose we keep the identification 9 — 9{R) in order to ensure that the zero temperature cancellation occurs. 
Then we see immediately that this cancellation cannot persist at finite temperature - this is because we always have 
the freedom to modify the profile function 9(r) in the external region, without changing its value at r = R. This will 
have no effect on the contribution from inside the bag because we have not changed 9 [which is equal to 9(R)], and 
we have not changed R. On the other hand, we see from (^) that the contribution from outside the bag does change. 
Therefore, the fractional parts of the internal and external contributions to the finite temperature fermion number 
cannot cancel. 

Thus, at finite temperature, the baryon number of the hybrid chiral bag is no longer unity for arbitrary configurations 
of the chiral field U. Physically, this says that the thermal expectation value (^) of the baryon number differs from 
its topological zero temperature value of 1 (which is a temperature independent vacuum polarization effect) due to 
the thermal occupation of excited quark states (which is a temperature dependent plasma response effect). The finite 
temperature corrections are nontopological because they are more sensitive to the details of the quark spectrum, 
which in turn are highly sensitive to the precise form of the chiral background. 

Finally, we mention that it would be interesting to investigate these issues analytically in the simpler cases of the 
(1+1) and (2+1) dimensional hybrid bags, where exact analytic information is available for the quark spectrum inside 
the bag [f^ |§. 
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